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SUMMARY 


The  study  of  the  supersonic  pressure  distribution  and  wave  drag  of 
sweptback  wings  at  zero  lift  has  been  extended  by  means  of  linearized 
theory  to  include  the  combined  effects  of  plan -form  taper  and  curvature 
of  profile.  -Calculations  have  been  made  for  a  specific  wing  with  a 
supersonic  leading  edge  and  these  results  are  applicable  to  a  limited 
family  of  plan  forms.  Comparisons  have  been  made  of  the  pressure  dis¬ 
tributions  and  the  wave  drag  calculated  by  this  exact  linear  solution 
with  the  results  of  an  arbitrary  simplified  approximation  based  on 
formulas  for  untapered  plan  forms  with  symmetrical  parabolic -arc  sections. 
In  addition,  a  comparison  of  the  wing  wave -drag  coefficient  has  been  made 
with  the  results  obtained  from  formulas  for  geometrically  similar  swept¬ 
back  tapered  wings  with  rhombic  sections  modified  by  a  profile -correction 
factor. 


INTRODUCTION 


During  the  past  few  years,  extensive  investigations  utilizing 
linearized  theory  have  been  in  progress  to  predict  the  aerodynamic 
characteristics  of  a  wide  variety  of  wing  configurations  at  supersonic 
speeds.  The  method  of  reference  1  was  applied  for  the  case  of  zero  lift 
in  references  2  and  3  to  predict  the  wave-drag  characteristics  of  4 
untapered  sweptback  wings  with  symmetrical  parabolic -arc  sections  and 
in  references  k  and  5  to  predict  the  wave  drag  of  tapered  sweptback  wings 
with  rhombic  profiles  (symmetrical  double-wedge  sections  with  maximum 
thickness  at  the  midchord).  Reference  6  considered  the  effect  of 
chordwise  location  of  maximum  thickness  on  the  supersonic  wave  drag  of 
sweptback  wings  with  symmetrical  double -wedge  profiles.  Since  many 
practical  configurations  include  combined  plan -form  taper  and  curvature 
of  profile,  it  is  of  current  interest  to  investigate  this  case.  This 
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paper  supplements  the  investigations  of  references  1  to  5  by  presenting 
the  calculations  for  the  pressure  distribution  for  a  limited  family  of 
tapered  sveptback  plan  forms  with  symmetrical  parabolic-arc  sections  and 
supersonic  leading  edges  at  zero  lift. 


SYMBOLS 


x,  y,  z 

V 


M 


Cartesian  coordinates  with  origin  located  at  wing  apex 
coordinates  of  origin  of  source  line 
free-stream  velocity  in  x-direction 
free -stream  Mach  number 

first  partial  derivative  of  disturbance -velocity 
potential  in  x-direction 


$xx>  0yy*  ^zz 


second  partial  derivatives  of  disturbance -velocity 
potential  in  x-,  y-,  and  z-directions,  respectively 


A 


b 


ct 

X 


*i>  yi 


c 

S 

A 

dz/dx 

m 


sweep  of  leading  edge 
wing  span 
root  chord 
tip  chord 

taper  ratio  (ct/cr) 

coordinates  of  intersection  of  extended  leading  and 
trailing  edges 


local  chord 


in  x-direction 


wing  area 

aspect  ratio  (b2/s) 
slope  of  airfoil  surface 

slope  of  source  line 
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Within  the  assumptions  of  small  disturbances,  the  linearized 
equation  for  the  velocity  potential  is 

(1  -  M2)0xx  +  0yy  +  0ZZ  =  0  (1) 

Because  of  the  linearity  of  equation  (l)  a  solution  may  be  used  to 
denote  one  of  the  velocity  components  rather  than  the  velocity  potential 
(See  reference  7-)  The  specification  of  one  component  in  this  manner 
actually  describes  the  whole  flow  field  since  the  other  components  may 
be  obtained  by  integrating  the  given  component  to  obtain  the  velocity 
potential  and  then  differentiating  the  results  along  the  desired 
directions  to  obtain  the  desired  components.  It  was  shown  in  refer¬ 
ence  1  that  a  solution  to  equation  (l)  corresponding  to  a  disturbance 
velocity  at  z  =  0  of  an  oblique  line  of  pressure  sources  in  the 
plane  z  =  0  is,  for  mp  >1.0,  equivalent  to 


b 
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and  for  m3  <  1.0, 

u 


where 


The  +  and  -  signs  in  equation  (4)  refer  to  equations  (2)  and  (3), 
respectively.  The  origin  of  the  lines  of  sources  described  in  equa¬ 
tions  (2)  and  (3)  is  the  origin  of  the  coordinate  system.  Equa¬ 
tions  (2)  to  (4)  were  given  in  references  2  and  3  for  a  line  of  sources 
with  origin  located  at  any  point  (  !,tj)  in  another  form  and  are  for  the 
supersonic  leading  edge  (m3  >  1) 


u6  (x,y)  =  R.P 


I  COS 


-1 


X  -  i  -  m3g(y  -  n) 

3  [m(x  -  I)  -  (y  -  t))_ 


and  for  subsonic  leading  edge  (m3  <  1) 


(5) 


ua,n(x'y) 


R.P. 


I  cosh-1 


x  -  I  -  m3g(y  -  t)) 
3|m(x  -  i)  -  (y  -  tj)| 


where  for  an  untapered  wing  m  =  cot  A  and  for  a  tapered  wing 


(6) 


m 


yi 

xi  -  s 


(7) 


In  reference  1,  the  following  is  stated:  "Curved  surface s  require 
a  continuous  distribution  of  sources  and  sinks  alined  with  the 
generators  of  the  surface.  Each  elementary  source  line  causes  an 
infinitesimal  change  in  direction  of  the  surface  and  hence  the  slope  at 
any  point  may  be  obtained  by  adding  up  the  effects  of  all  sources  ahead 
of  that  point."  Thus,  it  was  shown  (reference  1)  that 


dz  it  /  \/±( m232  -  1)  dl  t 

dx  v  Jx.  m  a  i a 


(8) 


and 


dl  =  V  m 

dE  "  1  l/±(m2p2  -  l) 


d2z 

d!2 


j m 


S 


* 


(9) 
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Pressure  Distribution 
The  pressure  coefficient  is  given  by 

P  =  -  (10) 

where  the  summation  in  equation  (10)  represents  the  superposition  of 
the  finite  sources  of  strength  sufficient  to  form  the  desired  angle  of 
intersection  of  the  arcs  of  the  leading  and  trailing  edges  together  with 
a  continuous  distribution  of  infinitesimal  sources  and  sinks  along  the 
chord  line  to  produce  any  desired  curvature. 

The  following  equations  from  references  3  and  2,  respectively, 
when  applied  to  wing  plan  forms  where  each  tip  affects  solely  its  own 
half  of  the  wing  can  be  given  as  follows:  For  a  supersonic  leading 
edge  (mQp  >  l) 


u  u0,0  +  u0,0  '  Du0,0  “Du0,0  *  ub/2mQ,-bf2  +  DUb/2mQ,b/2  ^ 
and  for  a  subsonic  leading  edge  (iiiqP  <  1 ) 

u  =  ^*0,0  +  *0,0  +  uc,0  +  *c,0  ”5U0,0  "  d*0,0  + 

Duc ,0  +  Duc,0  ~  ub/2mQ,b/2  +  Bub/2mQ,b/2  (12) 

where  the  unbarred  quantities  in  equations  (ll)  and  (12)  refer  to  con- 
tributions  of  the  sources  in  the  half- wing  where  the  point  under  con¬ 
sideration  is  located,  and  the  barred  quantities  refer  to  the  contri¬ 
butions  of  those  sources  in  the  other  half-wing,  the  disturbances  of 
which  are  propagated  to  that  point.  The  continuous  distributions  of 

sources  and  sinks  (represented  by  -ju  and  -gu)  in  equations  (ll)  and  (12 ) 

are  obtained  by  chordwise  integration  of  the  solution  for  an  elementary 
line  source  given  by  the  expression: 


-  B.P.  21  cos-1  iffa  -  ^  -  *7  -  g 


The  following  formulas  can  be  used  to  evaluate  the  contributions 
of  the  various  terms  of  equations  (ll)  and  (12):  For  m^p  >  1 


u0)0<x,y)  ♦  u0;0(x,y)  -  E.P.  21  cos’1 
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For  moP  <  1 


”0,0^)  +  uo^o^y) 


R.P. 


cosh' 


x  -  Jatfpy  ,__3  x  +  moP2y 

Pln^x  -  y  |  +  y] 


(15) 


For  the  region  ahead  of  the  Mach  line  from  the  wing  apex  (x  ^  py 
and  diqP  >  1) 


5uo,o(x,y)  +  liuo,o(x'y) 


(16) 


* 


The  contributions  ^~u  +  —uy  of  the  chordwise  distributions  of  source 

and  sink  lines  are  obtained  by  a  summation  of  the  effects  of  all  the 
elementary  source  and  sink  lines,  the  disturbances  of  which  are  propagated 
to  the  point  (x,y).  For  instance,  in  the  region  between  the  leading  edge 
and  the  Mach  line  from  the  wing  apex  (equation  (l6))  the  contributions 
of  the  leading -edge  source  lines  is  Ijt,  and  the  increment  due  to  the 
change  in  source  strength  from  an  infinitesimal  source  line  to  another 
approaches  jt  dl  in  the  limit.  In  the  limit,  the  summation  becomes 


an  integration  over  the  proper  limits  of 


*  dl 


or 


the  root  chord.  For  the  region  behind  the  Mach  line 
apex  (x  ^  Py  and  mQp  >  1) 


along 


from  the  wing 


D"o,o(x>y)’ +  D^o^y) 


+ 


£=x-py  d2z 

d£2 


cos 


-i  i  £ 


(x  -  S)2  -  pV 

(x  -  I)2  -  y2 


\A "  {§)' 


di 

(IT) 


» 
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where  m  is  a  function  of  £  given  by  equation  (7)*  The  first  term 
in  equation  (17)  is  the  contribution  of  the  elementary  source  lines 
from  the  leading  edge  to  the  Mach  line  emanating  from  the  wing  apex, 
and  the  last  term  represents  the  contribution  of  the  elementary  source 
lines  between  the  Mach  line  from  the  wing  apex  and  the  forward  Mach 
line  of  the  point  (x,y).  For  a  subsonic  leading  edge  ( mQp  <  l) 


r 


1^0,0  +  = 


M  &  C0sh-1  x  - «  -  -fr 

dg2 _ 3 im(x  -  g)  -  y 


1=0 


dl 


\ 


U 


a^z  co  . -i 


dl- 


l-0 


-  1 


'UlH±£  dS  l 


(18) 


J 


where  equations  (l6)  to  (18)  were  given  in  another  form  in  reference  1. 


In  calculations  for  a  parabolic -arc  section  in  the  stream  direc¬ 
tion,  d2z/d|2  is  constant  and  may  be  removed  from  the  integrands  of 
equations  (l6)  to  (18).  Because  of  the  complexity  of  the  integrations 
(equations  (l6)  to  (l8))  for  the  contribution  of  combined  plan-form 
taper  and  curvature  of  profile  to  the  pressure  distribution,  graphical 
integration  was  considered  the  most  advantageous  means  of  obtaining  the 
results. 


Drag 

Since  the  pressure  distribution  (equation  (10) )  for  the  tapered 
configuration  was  not  available  in  closed  form,  the  section  wave -drag 
coefficients  and  the  wing  wave -drag  coefficient  were  calculated  by 
graphical  means  from  the  equations: 

pi 


c  a  =  2 


o 


p  M 

dx  c 


(19) 
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where 


and 


M  =  1 £(k  .  *) 

dx  c\2  c  J 


(20) 


(21) 


RESULTS  AND  DISCUSSION 


Pres sure -distribution  calculations  have  been  made  for  the  wing 
configuration  of  figure  1  having  a  symmetrical  parabolic-arc  section. 

The  chordwise  pressure  distributions  at  the  root  chord,  tip  chord,  and 
four  intermediate  spanwise  stations  are  presented  in  figure  2.  The 
distribution  of  section  wave-drag  coefficient  (equation  (19))  is  shown 
in  figure  3  and  was  integrated  graphically  (equation  (21))  to  give  the 
wing  wave-drag  coefficient  Cp  =  5«3M t/c)^  for  the  wing  configuration 
of  figure  1  or  Q)  =  4.86(t/c)^cot  A  for  a  limited  family  of  tapered¬ 
wing  configurations  with  symmetrical  parabolic-arc  sections 
satisfying  X  =  0.531*  Ap  =  4.65*  ami  p  cot  A  =  1.375* 

As  a  point  of  interest,  the  pressure  distributions  were  calculated 
for  the  previously  mentioned  stations  by  an  arbitrary  approximation 
and  are  presented  in  figure  2;  the  pressure  distribution  at  each  station 
was  calculated  by  the  formulas  of  reference  1  for  untapered  wings  having 
the  same  span  and  leading-edge  sweep  as  the  tapered  wing  and  having  a 
constant  chord  equal  to  the  local  chord  of  the  tapered  wing  at  the 
spanwise  station  under  consideration.  The  untapered  wings  used  in 
these  computations  are  illustrated  in  figure  4.  Although  the  results 
compare  favorably  in  the  vicinity  of  the  root  chord,  the  agreement 
becomes  poor  with  distance  from  the  root  chord.  (See  fig.  2.) 

The  distribution  of  section  wave-drag  coefficient,  also  calculated 
with  the  use  of  the  geometric  assumptions  indicated  by  figure  4  and 
the  formulas  of  reference  3,  is  shown  in  figure  3.  The  corresponding 
wing  wave-drag  coefficient  is  Cp  =  5*29(t/c)^cot  A  *  which  is  about 
9  percent  higher  than  the  result  of  this  paper. 

The  wing  wave -drag  coefficient  calculated  in  the  present  paper 
(Cn  =  l(-.86(t/c)2cot  A)  was  compared  with  the  value  obtained  for  a 
geometrically  similar  wing  with  rhombic  profile  (reference  5)  multiplied 
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by  a  prof ile -correction  factor.  This  factor  was  chosen  as  the  ratio  of 
the  wing  wave -drag  coefficients  of  two  untapered  wings  (having  the  same 
aspect  ratio  and  leading-edge  sweep  as  the  tapered  wing)  with  parabolic- 
arc  (reference  3)  and  rhombic  sections  (reference  6),  respectively.  The 
wing  wave-drag  coefficient  Cp  =  4.85(t/c)^cot  A  so  obtained 
when  Ap  =  4.65  and  P  cot  A  =  1.375  compared  excellently  with  the 
value  calculated  herein. 


CONCLUDING  REMARKS 


The  study  of  the  supersonic  pressure  distribution  and  wave  drag  of 
sweptback  wings  at  zero  lift  has  been  extended  by  means  of  linearized 
theory  to  include  the  combined  effects  of  plan-form  taper  and  curvature 
of  profile.  Calculations  have  been  made  for  a  specific  wing  with  a 
supersonic  leading  edge  and  these  results  are  applicable  to  a  limited 
family  of  plan  forms.  Comparisons  have  been  made  of  the  pressure  dis¬ 
tributions  and  the  wave  drag  calculated  by  this  exact  linear  solution 
with  the  results  of  an  arbitrary  simplified  approximation  based  on 
formulas  for  untapered  plan  forms  with  symmetrical  parabolic -arc  sections. 
It  appears  that  these  approximations  are  adequate  only  in  the  vicinity 
of  the  root  chord.  In  addition,  the  wing  wave -drag  coefficient  was 
compared  with  a  value  obtained  from  that  of  a  geometrically  similar 
sweptback  tapered  plan  form  with  a  rhombic  profile  by  using  as  a  profile - 
correction  factor  the  ratio  of  the  wing  wave -drag  coefficients  of  two 
untapered  wings  (having  the  same  aspect  ratio  and  leading-edge  sweep  as 
the  tapered  wing)  with  parabolic -arc . and  rhombic  sections,  respectively. 
This  approximation  indicated  excellent  agreement  for  the  cases  considered. 
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Figure  2.-  Illustrative  chordwise  pressure  distribution  on  a  symmetrical 
parabolic-arc  section  at  zero  lift,  mg 3  =  1.375. 


Untapered-wing  formulas  (reference  3) 


Figure  3--  Wave-drag  distribution  on  a  symmetrical  parabolic-arc  section  at  zero  lift,  mQp 


Figure  4.-  Method  of  replacing  sections  of  tapered  wing  by  corresponding  sections  of  untapered  wings. 
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